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We use ah initio simulations to investigate the incorporation of pyridine molecules (C5H5N) in the van 
der Waals gaps of Bi2Se3. The intercalated pyridine molecules increase the separation distance between 
the Bi2Se3 quintuple layers (QLs), suppressing the parity inversion of the electronic states at the T-point. 
We find that the intercalated region becomes a trivial insulator. By combining the pristine Bi2Se3 region 
with the one intercalated by the molecules, we have a non-trivial/trivial heterojunction characterized by the 
presence of (topologically protected) metallic states at the interfacial region. Next we apply an external 
compressive pressure to the system, and the results are (i) a decrease on the separation distance between the 
QLs intercalated by pyridine molecules, and (ii) the metallic states are shifted toward the bulk region, turning 
the system back to insulator. That is, through a suitable tuning of the external pressure in Bi2Se3, intercalated 
by pyridine molecules, we can control its topological properties; turning-on and -off the topologically protected 
metallic states lying at the non-trivial/trivial interface. 


Three dimensional topological insulators (TIs) are in¬ 
sulator materials in the bulk phase, but they present 
metallic topological surface states (TSSs), which are pro¬ 
tected by time-reversal symmetry; as a result backscat- 
tering processes by time-reversal invariant impurities or 
defects are forbidden^. These materials are of great 
promise for spintronics applications, due to the formation 
of nearly dissipationless spin-polarized surface current^. 
Currently, Bi2Se3 is one of the most investigated TI due 
to its large band gap (^ 0.3 eV), nearly idealized single 
Dirac cone and for being a pure compound^”^. Bi2Se3 
presents a rhombohedral structure composed by quintu¬ 
ple layers (QLs) of Se and Bi atoms, forming a sequence 
of Se-Bi-Se-Bi-Se atoms covalently bonded; these QLs 
are stacked along the c-axis of a hexagonal structure by 
van der Waals (vdW) interactions. 

Recent experiments have been exploring the possibil¬ 
ity of inserting guest species in the vdW gaps of Bi2Se3, 
this process is known as intercalation. Koski et al. 
have added various zerovalent metals in the vdW gaps of 
Bi2Se3 nanoribbons, it is expected that new properties 
and/or tuning of the Bi2Se3 properties can be achieved 
by intercalation, e.g. superconductivity^. Beside met¬ 
als, molecules have been intercalated in Bi2Se3, for in¬ 
stance, it has been shown that pyridine molecules in 
Bi2Se3 present interesting properties for optoelectronic 
applications ^. 

In this work we aim to investigate the geometry and 
electronic structure of Bi2Se3 intercalated by pyridine 
molecules (py-Bi2Se3). Upon the presence of the pyri¬ 
dine molecules between the QLs, we find that the py- 
Bi2Se3 system becomes a trivial insulator. By consider¬ 
ing a hetero junction composed by py-Bi2Se3 and pristine 
Bi2Se3 (py-Bi2Se3/Bi2Se3), we analyze (i) the occurrence 
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of topologically protected metallic states embedded in 
py-Bi2Se3/Bi2Se3, and (ii) the control of those metallic 
states (turning-on and -off) upon the application of ex¬ 
ternal pressure. 

The calculations are performed based on the density- 
functional theory (DFT) as implemented in the Vienna 
ab initio simulation package (VASP)^^. We use the gen¬ 
eralized gradient approximation (GGA), in the form pro¬ 
posed by Perdew, Burke and Ernzerhof^^, to describe 
the exchange-correlation functional. The Kohn-Sham or¬ 
bitals are expanded in a plane wave basis set with an en¬ 
ergy cutoff of 400 eV. The electron-ion interactions are 
taken into account using the Projector Augmented Wave 
(PAW) method^^. All geometries have been relaxed un¬ 
til atomic forces were lower than 0.02 eV/A. The Bril- 
louin Zone is sampled according to the Monkhorst-Pack 
method^^, using at least a 3x3x1 mesh. We have also 
used a functional that accounts for dispersion effects, rep¬ 
resenting van der Waals (vdW) forces, according to the 
method developed by Tkatchenko-Scheffler (TS)^^, which 
is implemented on VASP^^. The inclusion of van der 
Waals forces in the simulations is necessary to obtain the 
correct vdW gap between two QLs^^, the interaction be¬ 
tween molecule and Bi2Se3 is also better described with 
the inclusion of vdW interactions. 

Initially we calculate the equilibrium distance (zq) be¬ 
tween two consecutive QLs of Bi2Se3, vdW gap. Here we 
have considered two isolated QLs and minimize the total 
energy as a function of the width of vdW gap (Az) 

As shown in Fig. 1 (squares), we find the energy mini¬ 
mum for zo = 2.65 A; this result is in agreement with ex¬ 
perimental measurements and recent theoretical stud¬ 
ies I6,19 gy using the same procedure, we determine zq 
upon the presence of a pyridine molecule intercalated in 
between the QLs. The QLs present a 2x2 surface period¬ 
icity, which corresponds to a molecular concentration of 
1.64 xlO ^ mol./A^. In this case, we find that (i) there 
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FIG. 1. Formation energy for two Bi 2 Se 3 QLs at various 
separation distances. Black line (circles) represents the pris¬ 
tine system and red line (squares) the system with a pyridine 
molecule intercalated between the 2 QLs. 


is no chemical bonding between the molecule and the 
Bi2Se3 QLs, and (ii) the vdW gap increases by 3.85 A, zq 
= 2 . 65 ^ 6.50 A, Fig. 1 (circles). Such an increase of the 
vdW gap reduces the binding energy between the QLs, 
however, the formation of pyridine intercalated Bi2Se3 
QLs is still an exothermic process. 

TIs present a parity inversion between valence and 
conduction bands due to spin-orbit coupling (SOC). For 
Bi2Se3 the bulk phase presents a band parity inversion 
at the F-point, namely the SOC induces a band inver¬ 
sion between Se p2;-orbitals in the valence band and Bi 
p^-orbitals in the conduction band^. We use this charac¬ 
teristic to determine whether the material is behaving as 
trivial or topological insulator. Thus, by turning-off the 
SOC, we observe that higher energy levels of the valence 
band have higher contribution from the Se p2;-orbitals, 
while the lower energy levels of the conduction band are 
dominated by Bi p2;-orbitals. Meanwhile, by turning-on 
the SOC we observe an inversion between valence and 
conduction band orbitals around the F-point, promot¬ 
ing the (band) parity inversion. The same procedure was 
applied to examine if such a band inversion still occurs 
in py-Bi2Se3. By considering the pyridine concentration 
of 1.64 xlO ^ mol./A^ intercalated in Bi2Se3, we verify 
that when the SOC is not included the py-Bi2Se3 sys¬ 
tem is an insulator, where the valence band maximum 
is dominated by Se p2;-orbitals and the conduction band 
minimum by Bi p2;-orbitals [Fig. 2 (a)], as observed for the 
pristine Bi2Se3 bulk. In contrast, by turning SOC on we 
observe the lack of band inversion around the F-point 
[Fig. 2 (b)], suggesting that pyridine molecules interca¬ 
lated in Bi2Se3 bulk suppress its topological properties, 
namely, the py-Bi2Se3 system is a trivial insulator. 

This topological-trivial transition in the Bi2Se3 struc¬ 
ture is possibly caused by the increase in the vdW gaps 
due to the pyridine intercalation. In order to verify such 
an assumption, we calculate the electronic band struc¬ 
ture of Bi2Se3, but keeping the equilibrim geometry of 



FIG. 2. Band structures for (a) Bi 2 Se 3 bulk with intercalated 
pyridine molecules excluding SOG, and in (b) including SOG. 
(c) Bi 2 Se 3 bulk with the vdW gap obtained in the system with 
intercalated molecules (Fig. 1), excluding SOG; and in (d) 
including SOG. The symbols are contributions of pz-orbitals 
for Se (red) and Bi (green) atoms. 


the pyridine intercalated Bi2Se3 system. Our results are 
presented in Figs. 2 (c) and 2 (d), by turning-off and -on 
the SOC we observe that the band inversion at the F- 
point is still absent. 

In experiments we would most likely observe regions 
with pyridine molecules intercalated between the vdW 
gaps and regions of pristine Bi2Se3, as a result we 
may have a trivial/topological insulator heterojunction 
composed by py-Bi2Se3 and Bi2Se3 (py-Bi2Se3/Bi2Se3). 
To simulate such a hetero junction, we have consid¬ 
ered a supercell containing 9 QLs, where we have three 
QLs intercalated by the molecules and six consecu¬ 
tive pristine QLs, (py-Bi2Se3)3/(Bi2Se3)6. As shown in 
Fig. 3 , the pyridine molecules are intercalated between 
QL(1) and QL( 4 ), characterizing the (py-Bi2Se3)3 re¬ 
gion, and QL( 4 )/QL( 5 )/.../QL( 9 )/QL(l) forms the pris¬ 
tine (Bi2Se3)6 region. Given that geometry, the triv¬ 
ial/topological interface are characterized by QL(1) and 
QL( 4 ). Band structure calculations reveal the presence 
of metallic states forming a Dirac-like cone at the F- 
point, indicated by solid black lines in Fig. 3 . To lo¬ 
calize these states we compute the pz-orbitals contribu¬ 
tion from various QLs around the F-point, and project 
these states in the energy levels shown in the band struc¬ 
ture^^. We find that the electronic states around the 
Dirac point are mainly attributed to the Bi p^-orbitals 
lying at the py-Bi2Se3/Bi2Se3 interface region, QL(1) 
and QL( 4 ). In contrast, there are no electronic contri¬ 
bution to the metallic states from the py-Bi2Se3 bulk 
region, QL(2) and QL( 3 ). Somewhat similarly, the elec¬ 
tronic contributions to the Dirac-like cone coming from 
the bulk region of Bi2Se3 [QL( 5 )-QL( 8 )] are almost neg¬ 
ligible. Those findings allow us to conclude that (py- 
Bi2Se3)ni/(Bi2Se3)n heterostructures (or superlattices) 
give rise to topologically protected metallic states, em- 
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FIG. 3. Each panel shows the band structure (full black lines) 
for the structure shown on the right side. The pz-orbital 
states is projected on the band structure for the first four QLs, 
where green circles represent Bi pz-orbitals and red circles Se 
pz-orbitals. 


bedded at the trivial/topological interface region. How¬ 
ever, it is worth noting that the appearence of those topo¬ 
logically protected sufaces states depends on the topo¬ 
logical film (Bi2Se3) thickness (n). In Ref. the au¬ 
thors verified the presence of converged TSS for 3 QLs 
of Bi2Se3. On the other hand, in a recent study, we ver¬ 
ify the formation of TSSs for an interlayer spacing, i. 
e. vdW gap, larger than ^ 5.5 A Thus, we can infer 
that it is possible to get topologically protected metallic 
channels by considering just a single layer of py-Bi2Se3 
embedded in Bi2Se3. 

In order to verify the statement above, we construct a 
supercell containing one molecule separated by six Bi2Se6 
QLs, (py-Bi2Se3)i/(Bi2Se6)5, this structure is shown in 
Fig. 4 (c). Indeed, we find the TSSs forming a Dirac- 
like cone at the F-point, as shown in Fig. 4 (a). Those 
states lie at the (py-Bi2Se3)i/(Bi2Se3)5 interface region. 
The spin-texture of TSSs is constrained by the time re¬ 
versal symmetry. Here we calculate the expected value 
of spin-polarization components ((5'^(k))) for the TSSs, 
of (py-Bi2Se3)i/(Bi2Se3)5, near the Dirac point. (5'n(k)) 
can be written as = (^/2)(0n(k)|crc^|0n(k)), in 

cartesian coordinates {a = x^y^z), where a represents 



FIG. 4. (a) Band structure for the system depicted in (c), 
along the path shown in the projected surface 2D Brillouin 
zone represented in (b). The inset in (a) shows the spin- 
texture inside the blue square region around the F-point. 


the Pauli matrices, and (/)n(k) the single particle Kohn- 
Sham eigenfunction, and n represents the band index. 
We have considered 0n(k) for wave vectors (k) along 
the F-M and F-K directions. As depicted in the inset of 
Fig. 4 (a), for the electronic states parallel to the F-M di¬ 
rection, we find (i) positive (negative) values of ( 5 'n,cc(k)) 
for the occupied (empty) states, while (ii) {Sn,y(k)) = 
(5'n^;2(k)) = 0 ; whereas for k parallel to the F-K direc¬ 
tion, we find (iii) negative (positive) values of {Sn,y{^)) 
for the occupied (empty) states, (iv) ( 5 'n,cc(k)) = 0, and 
(v) {Sn,zO^)) values are negligible around the F-point. 
Such a picture of {Sn,a 0 ^))^ for the TSSs lying at the 
(py-Bi2Se3)i/(Bi2Se3)5 interface, is practically the same 
as that obtained for the TSSs on the Bi2Se3(lll) sur¬ 
face 

We now propose a mechanism to gain control over 
the topologically protected metallic states in (py- 
Bi2Se3)m/(Bi2Se3)n. In particular, we considered the 
(py-Bi2Se3)i/(Bi2Se3)5 system. We apply a compres¬ 
sive pressure directed along the c-axis, namely normal 
to the QLs, decreasing the separation distance in the 
vdW gap (Az) from 6.50 A to 4.12 A, upon an external 
pressure of P « 9.1 GPa. This (pressure) value is com¬ 
parable to the ones achieved experimentally for Bi2Se3 
structures^^’^^. The energy difference between the re¬ 
laxed and compressed (py-Bi2Se3)i/(Bi2Se3)5 is shown 
in the upper part of Fig. 5 . An energy of ^ 10 meV/A^ 
needs to be added to compress the system from Az = 
6.50 to 4.12 A. In the lower part of Fig. 5 we show the 
band structure for four values of Az, starting from 6.50 A 
which has already been shown to have metallic states. 
For Az = 5.65 A, the metallic states at the F-point are 
still present. Upon further decrease of A^, the metallic 
states start to move in the bulk band direction. For Az 
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FIG. 5. Upper panel: Energy per area to compress the re¬ 
laxed structure shown in Fig. 4 (c), reducing the vdW gap 
between the QLs separated by the molecule. Lower panels: 
Band structure for the various Az values. 


= 4.71 A, the Dirac point has been suppressed, as we 
find an energy gap at the F-point. The energy gap at the 
F-point is even larger for Az = 4.12 A, washing out the 
TSSs at the (py-Bi 2 Se 3 )i/(Bi 2803)5 interface region. By 
removing the pressure the system can reversibly return to 
its original geometry, and again the metallic states will be 
present. Thus, we can combine pressure application and 
pyridine molecules intercalation to the Bi 2 Se 3 structure 
to turn the topologically protected metallic states on and 
off, going from an insulator to a semi-metallic material 
by taking advantage of the TI properties of Bi 2 Se 3 . The 
results presented above can be extended to any number 
of QLs intercalated by molecules, assuring that we have 
a Bi 2 Se 3 /py-Bi 2 Se 3 junction. 

We have performed an ab initio study of pyridine 
molecules intercalated in the vdW gaps of Bi 2 Se 3 (py- 
Bi 2 Se 3 ). In py-Bi 2 Se 3 the inter-QL distance increases, 
which turns the TI material into a trivial insulator. By 
considering (py-Bi 2 Se 3 )ni/(Bi 2 Se 3 )n heterojunctions, we 
find a trivial/topological interface, characterized by the 
presence of topologically protected metallic states (form¬ 
ing a Dirac-like cone) embedded at the interface region of 
the heterostructure. Such metallic states can be present 
even for a single QL incorporated by pyridine molecules, 
(py-Bi 2 Se 3 )i/(Bi 2 Se 3 )n. Lastly, we have shown the pos¬ 
sibility to control the occurrence of such metallic states 
in (py-Bi 2 Se 3 )m/(Bi 2 Se 3 )n, upon an external compres¬ 
sive strain; turning-on and -off those metallic states at 
the hetero junction interface. 
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